OPEN HEAVY FLAVOR PRODUCTION 

IN DEEPLY INELASTIC ep SCATTERING AT HERA* 

KARIN DAUMt 
Universitdt Wuppertal, Rechenzentrum, Gaufistrafie 20 
D-42097 Wuppertal, Germany 
E-mail: daum@mail.desy.de 

(On behalf of the HI and ZEUS collaborations) 

"t> ' 

ON . ABSTRACT 
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Recent results on inclusive D° and D*^ production in deeply inelastic ep 
bJQ' scattering at y/s = 301 GeV are summarized. The data have been collected 

by the HI and ZEUS experiments at HERA. The total and the differential 
cross sections are discussed in the framework of LO and NLO QCD predictions. 
| The data exhibit clear evidence for boson gluon fusion dominating open heavy 

flavor production in the kinematic range currently explored at HERA. The 
t— I | measurements of Fg c (x,Q 2 ) at small Bjorken x are presented. The prospects 

for future analyses of open charm and beauty production including detector 
upgrades and anticipated high luminosities are investigated. 

o : 

1. Introduction 

o 

■ 

Q\ • The study of heavy flavor production in deeply inelastic lepton nucleon scattering 

(DIS) is of importance for the understanding of the parton densities in the nucleon 0BIB 
In accordance with early measurements by the EMC collaboration i, neutral current 
charm quark production is expected to be dominated by the boson gluon fusion (BGF) 
process ^f{Z Q )g — > cc at large Bjorken x. Recent results on inclusive production of 
f> ' D° and D*^ mesons from HI I and ZEUS i at HERA enable the charm production 

\ mechanism at small x to be investigated. 

The measurement of the charm contribution F^(x, Q 2 ) to the proton structure 
function F2 is supposed to provide information on this mechanism. Present parton 
density calculations include heavy flavors in the proton using different assumptions. 
Charm is either considered as flavor excitation (FE) ,i.e. ■j(Z°)q — > qX, in the 
massless quark evolution starting with a charm density of zero at a scale Qa ^ 
or it is exclusively treated as BGF taking into account the charmed quark mass tL2l. In 
terms of phenomenological inclusive analyses 100 BGF is predicted to dominate near 
threshold, while FE should lead to a good description well above threshold. Various 
schemes to match these two regions have been proposedi'0. Both approaches lead to 
similar predictions for F2 C (x, Q 2 ) but differences in exclusive heavy quark spectra are 
expected. 

If charmed hadron production is dominated by BGF, it provides direct information 
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Fig. 1. (a) The K 7r + mass distribution observed in the D analysis of DIS events from HI, and 
(b) the Am distribution obtained in the D* + analysis of DIS events from ZEUS. 



on the gluonic content of the proton, x g g(x g , Q 2 ), complementary to what is obtained 
from scaling violation analyses of F 2 . Since the reconstruction of hadrons with heavy 
flavor content unequivocally tags the parent quark, heavy quark spectra become 
directly accessible. Compared to two jet studies the analysis of charmed hadrons 
allows (a) to investigate the BGF process down to the phase space limit s = Am 2 c 
which extends the range for the determination of the gluon density in the proton 
significantly towards smaller x g , and (b) to tag the gluon in the proton almost free 
of background from FE of light quarks in the proton, since heavy flavor production 
from gluon splitting is suppressed by heavy quark masses. 

2. Results 

Published results on charm production in deeply inelastic ep scattering exist from 
both experiments HI and ZEUS based on an integrated luminosity of roughly 3 pb _1 
collected with each experiment at HERA in 1994. In addition preliminary results from 
ZEUS are available based on an integrated luminosity of 6.4 pb _1 from the HERA 
1995 running. The HI collaboration I has performed the tagging of heavy quark events 
by reconstructing Z)°(1864) a and D* + (2010) mesons, while the ZEUS collaboration 
gives results for the inclusive D* + (2010) analysis. 

The D° is identified via its decay mode D° —>■ K~ir + and the D* + through the 
decay chain D* + — > D°7i^ low — > K~7i + 7i^ low . For the latter use is made of the low 
Q-value in the decay of D* + — > D nf low which leads to a better resolution for the 
mass difference Am = m(D n^ low ) — m(D°) than for the D* + mass itself. 

Figure |I| shows the m^n distribution obtained in the inclusive D° analysis of HI 
and the Am distribution as observed in the inclusive D* + of ZEUS. Evidently the 
number of observed events containing heavy quarks is small. Only of the order of 100 
to 200 charmed mesons are identified in any of the different analyses. 

"Charge conjugate states are always included. 
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Fig. 2. Normalized xd distribution in deep inelastic ep scattering at (W) ~ 125 GeV for \t]d\ < 1-5. 
The open/closed points represent the D°/D* data of HI. See text for further explanations. 



2.1. Mechanism of charm production in DIS 

Due to the hard fragmentation of D mesons its momentum is strongly correlated 
with the parent charm quark momentum. Therefore the distribution of the scaled D 
momentum in the j*p system Xd = 2|Pp|/W depends on the charm quark momentum 
spectrum and thereby on the mechanism of charm production. In BGF a cc pair is 
produced recoiling against the proton remnant in the j*p system, whilst in FE a 
single charmed quark is emitted opposite to the proton remnant, confining the second 
charmed quark necessary for quantum number conservation. Consequently the mean 
(xd) of D mesons in the current and central fragmentation region is expected to be 
roughly a factor of 2 smaller for BGF compared to FE. 

Figure ||] shows the distributions \jo der/dxp obtained from the D° and D* + 
analysis of HI in comparison to the expectation for the BGF and FE process. The 
LO BGF prediction from AROMA H (shaded histo gram) agrees well with the data. 
The figure also includes the FE expectations for charmed mesons either by using the 
LEPTO 6.1 generator 0, from which only FE events are selected (dashed histogram), 
or by extrapolating the results from charm production in charged current v N 
scattering EH to HERA energies (full line). Large differences in the shape are observed 
for the data and these FE expectations. From a fit to the data of BGF and FE in 
arbitrary proportions it is concluded that more than 95% of neutral current charm 
production in DIS is due to the BGF process at (Q 2 ) ~ 25 GeV. This observation 
seems to contradict recent inclusive calculations of charm production in DIS §, from 
which it was concluded that for the kinematic range of the current analyses at HERA 
charmed quarks may already be considered as massless partons in the proton. 

2.2. Differential inclusive D* cross sections 

Figure || shows the differential e + p — > D*^X cross sections as a function of a) 
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b) ii(D*), c) W and 

d) Q 2 from ZEUS. The inner/outer error bars correspond to the statistical/total errors respectively. 
The band and histogram are explained in the text. 



Pt(D*), b) r](D*), c) W and d) Q 2 measured by ZEUS. The shaded band represents 
the result of the analytic NLO calculation! using the GRV-HO 94 gluon density of 
the proton which assumes charm being exclusively produced by BGF and treats charm 
quarks as massive particles. Both the renormalization scale and the factorization scale 
are set to /x = J Q 2 + 4m^. The largest uncertainties in the calculations arise from 
the ignorance of m c . The upper (lower) limit corresponds to a charm quark mass 
m c = 1.35(1.7) GeV. This calculation reproduces well the shapes of the differential 
distributions. The figure also includes the analytic LO calculation as histograms. The 
LO calculation describes the shapes of the distributions equally well. 

2.3. Integrated charm cross sections 

All differential distributions discussed in the previous sections are found to be 
well described by the BGF pro cess. Therefore, it is used to extrapolate the inclusive 
D cross sections, a(ep — > eDX), from the visible kinematic region to the full phase 
space in order to determine the charm production cross section, a(ep — > eccX), in 
DIS. 



Table 1. Measured and predicted integrated charm cross sections a(ep — > eccX) for different bins 
in Q 2 in the region y < 0.7 (ZEUS) and 0.01 < y < 0.7 (HI). 





m c [GeV] 


5 < Q 2 < 10 GeV 2 


10 < 


Q 2 < 100 GeV 2 


ZEUS D* ± 




13.5 ±5.2 ±1.8+};^ nb 


12.5 


±3.1 ±1.8+};^ nb 


LO MC 


1.5 


12.5 nb 




14.2 nb 


LO GRV 


1.5 


11.0 nb 




12.4 nb 


NLO GRV 


1.5 


9.4 nb 




11.1 nb 


HI D*± 






15.1 


± 1.8 + 2o ± 1.2 nb 


HI D° 






19 


5±2.6l 2 : ™nb 


NLO MRSH 


1.3 






11.3 nb 


NLO MRSH 


1.7 






8.1 nb 


HI F 2 fit 


1.5 






13.6 nb 



Table [I] summarizes the experimental results together with the predictions of 
the NLO calculations ii for different parton densities and different values of m c . 
The extrapolation in phase space is calculated in NLO using the GRV-H0 gluon 
density in case of the ZEUS data while the acceptance for the HI data is calculated 



using the MRSH gluon density. In both cases m c = 1.5 GeV and /i = \JQ 2 ± 4m 2 
are chosen. The charm fragmentation is described by the Peterson function with 
e c = 0.035±0.009 (ZEUS) and the symmetric Lund function with the standard setting 
of JETSET 7.4 (HI). The errors correspond to the statistical error, the experimental 
systematic uncertainty and the model dependent uncertainty respectively. The latter 
accounts for the possible variation of the relevant parameters of the model used for 
extrapolation, i.e. xg(x,Q 2 ), m c , /i and the fragmentation function. In the region of 
overlap the cross section inferred from the D*^ analyses are in reasonable agreement. 
The cross section from the D° analysis is higher but the few statistics do not allow 
any conclusion to be drawn. 

In table |l] the data is also compared with predictions from analytic NLO calculations 
for different gluon densities and different values of m c . Also given are LO analytic 
calculations and the result from the LO Monte Carlo. In addition HI has performed 
a NLO calculation using the gluon density as derived from the analysis of scaling 
violations of F 2 0. Within reasonable variations of the most relevant parameter, m c , 
the NLO calculations tend to be below the experimental results except for the gluon 
density obtained from the NLO F 2 fit. The LO analytic calculation and the LO MC 
yield cross sections generally in better agreement with the data. 
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Fig. 4. TTie charm contribution F^ c to the proton structure function as derived (a) from the D* + 
(dots) and D° analysis (squares) by the HI collaboration and (b) from the D* + by the ZEUS 
collaboration from the 1994 data (dots) and the 1995 data (triangles). The inner/outer error bars 
refer to the statistical/total errors. The EMC data are also shown (open boxes). The shaded band 
represent the predictions from the NLO fits to the HI and ZEUS F2 measurements respectively. The 
ZEUS 1995 data and the prediction of the NLO fit to the F2 measurement are preliminary. In (a) 
predictions from NLO calculations based on GRV-HO (full line), MRSH (dashed line) and CCFM 
(dotted line) gluon densities using a charmed quark mass of m c = 1.5 GeV are included 



2.4- Charm contribution _F 2 CC to the proton structure 

The charm contribution F| c to the structure function is obtained by using the 
expression for the one photon exchange cross section for charm production 

The contribution from Z exchange is expected to be small in the kinematic range 
explored presently. The _F£ C contribution has been estimated to be negligible 0J— 3. 
Therefore both F^ c and F£ c are neglected. 

The F£ c measurements from HI and ZEUS are displayed in figure [| together 
with the results from the EMC collaboration at large x. Within the limited statistics 
the results on F 2 CC from HI and ZEUS are in good agreement. The measurements at 
HERA extend the range of the F£ c measurements by two orders of magnitude towards 
smaller x values. The comparison of the HERA and EMC data reveals a steep rise of 
F 2 CC with decreasing x. HI has derived a ratio of F^jFi of 

Ff/F 2 ) = 0.237 ± 0.021 ± 0.041 , (2) 



which is an increase of about one order of magnitude of the overall charm contribution 
compared to the EMC result. This is consistent with the measured rise of the gluon 
distribution towards low x and the dominance of the BGF process observed here. 
Within the accuracy of the charm data this ratio seems to be constant in the kinematic 
range explored in the HI analysis. 

The data are also compared to the NLO calculations of F 2 C based on the gluon 
densities extracted from the NLO fits to the inclusive F 2 . In figure |]a the prediction 
from the HI F 2 data, using a charmed quark mass of m c = 1.5 GeV, is shown by 
an error band which includes the propagation of the statistical and the uncorrelated 
systematic errors on the total F 2 data through the fitting procedure. The uncertainty 
in the prediction of F 2 CC due to variations of m c is indicated by the error band in figure 
(|b for the ZEUS F 2 data. The upper and lower limits correspond to m c = 1.35 GeV 
and m c = 1.7 GeV respectively. No error propagation in the fit is performed. In the 
region of overlap in Q 2 the HI result of this fit favors a slightly larger gluon density 
at small x compared to the ZEUS measurement. 

Figure |]a also includes the predictions from the NLO calculation based on GRV- 
HO, MRSH and CCFM 13 gluon density calculations. Although the procedures for 
treating heavy flavors and for determination of the the gluonic content in the proton 
are quite different for these three parton density calculations, the resulting predictions 
for F| c are very similar. The accuracy of the data in terms of statistical and systematic 
uncertainty is by far not sufficient to be able to get more insight on the gluon in the 
proton from purely inclusive charm measurements. 

3. Future 

Although the charm production cross section is found to be about 25% of the 
total DIS cross section at HERA, the number of identified charm events is fairly 
small at present. For the results presented here heavy flavor tagging is performed 
by the reconstruction of charmed hadrons. From the experimental point of view one 
advantage of this method is the small contribution to the systematic uncertainties due 
to hadronization effects. Furthermore, as a consequence of the hard fragmentation of 
charm quarks, the reconstructed charmed hadron gives direct access to the parent 
charm quark. This allows the study of the production dynamic of heavy quarks almost 
free of background from other processes. The drawback of this tagging method is the 
small tagging efficiency (e tag ~ 0(1CT 3 )) due to the small probability that a heavy 
quark fragments into a specific heavy flavored hadron which subsequently decays into 
the desired decay mode, e.g. for the chain c — > D X — > D°irfX — > (K T 7T ± )-KfX 
this probability is only 6.7 ■ 10~ 3 . 

In order to increase et ag of heavy flavor production a Silicon Vertex Detector 
is installed in HI this year (1997) and will be installed in ZEUS together with 
the luminosity upgrade of the HERA machine in 1999/2000. Such a device enables 
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Fig. 5. The charm contribution FSf as expected for an integrated luminosity of Lint ~ 200 pb 1 using 
fully reconstructed charmed hadrons together with the Silicon Vertex Detector. 

the reconstruction of the decay vertices of heavy flavored hadrons and will thereby 
suppress the combinatorial background in the mass reconstruction of heavy flavored 
hadrons. The reduction of background will give access to a large variety of decay 
modes with good signal to background ratios. Considering the analysis of all accessible 
decay modes together e tag will be increased by about one order of magnitude. 

Alternatively heavy flavor events may be tagged by secondary vertices (SVTX) 
having high significance without explicit reconstruction of heavy flavored hadrons. 
This method will lead to tagging efficiencies comparable or even superior to the 
efficiency expected from the mass reconstruction method. Due to the absence of 
a mass constraint the SVTX method is expected to have larger background from 
faked secondary vertices. A much better understanding of fragmentation, detector and 
reconstruction performance is needed in this case to achieve experimental systematic 
uncertainties comparable to those of the mass reconstruction method. Furthermore it 
will become more difficult to access exclusive quantities of the production dynamics 
since the strong correlation between the heavy flavored hadron and the parent quark 
gets lost or will be weakened in partly reconstructed heavy flavored hadrons. 

Figure |5| shows the expected accuracy of _F 2 CC measurements for an integrated 
luminosity of Ci n t ~ 200 pb -1 using fully reconstructed charmed hadrons together 
with the Silicon Vertex Detector. The data points show the expectation based on the 
gluon density extracted from the HI NLO fit to the total F 2 for m c = 1.5 GeV and 
yU 2 = Q 2 + Am 2 c . The inner/outer error bars refer to the statistical/total errors. Also 
shown are the NLO predictions for the GRV-HO (full line) and MRSH (dashed line) 
parton densities and the same charm mass and scales. At low Q 2 the measurement 
will be limited by the experimental uncertainty while at large Q 2 the precision will be 
limited by statistics. The ultimate precision of 10 % will allow these parameterizations 
of the gluon density to be distinguished. 
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Fig. 6. The ratio F^/F^ as expected for an integrated luminosity of Ci n t ~ 500 pb 1 using fully 
reconstructed charmed hadrons together with the Silicon Vertex Detector. 



Figure ^ shows the expected accuracy of the ratio F^/F^ for an integrated 
luminosity of Ci nt ~ 500 pb -1 using fully reconstructed charmed hadrons together 
with the Silicon Vertex Detector @. Even for a luminosity of Ci n t ~ 500 pb" 1 the 
measurement of this ratio will be limited by statistics in the full x and Q 2 plane. 
For a given Q 2 or x the ratio is expected to rise with decreasing x or increasing Q 2 , 
respectively. A mean value of (F 2 66 /F 2 CC ) w 0.02 is predicted. 



4. Conclusions 



First results on charm production in deeply inelastic ep scattering at HERA have 
been presented. The investigation of differential D cross sections has shown that the 
production dynamics of charmed mesons in the current and central fragmentation 
region may be described by the boson gluon fusion process. The contribution from 
flavor excitation of charm in this region is found to be negligible at (Q 2 ) ~ 25 GeV 2 . 
This observation sets constraints on the kinematic region where charm may be treated 
as massless parton in the proton. 

The measurement of FSfix, Q 2 ) at small x reveals a strong rise of F^ix^ Q 2 ) from 
the fixed target to the HERA regime. Agreement is observed between the measured 
F 2 CC and the results from the NLO fits to the inclusive F 2 . Averaged over the kinematic 
range currently explored a ratio ^F 2 CC /F 2 ^ = 0.237 is observed which is one order of 
magnitude larger than at large x. 

An estimation on the accuracy for the exploration of charm and beauty contribution 
F| c and F^ b to the structure of the proton for an integrated luminosity of more than 
200 pb -1 has been summarized. For Q 2 < 100 GeV 2 the charm analysis will be 
systematically limited. The heavy flavor data will allow a direct measurement of the 
gluon density in the proton to be performed. A precision on xg(x, Q 2 ) of 10 % may 



be achieved from charm data, which is competitive with other methods. Despite of 
the smallness of the ratio F$ b / F^ c high luminosity running will allow detailed studies 
of beauty production dynamics. 
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